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The unstable burning of a biprovellent rocket comiustion 
chenber is investiested end a study meade of the resuiromonts for an 
entomatic closed loot control cirenit to stabilize the motor. 

The biprovellant combustion chomber ecrations dovelopod by 
Dre Le Croce 2} are utilised as the analytical descrirtion of the 
rocket motor burning phenomena. Honations similar to those devaloned 
by Dr. H. 8. teten'?? are used for the oxidizer and Mel surply systems 
and the tyo closed loop stebilizine circuits. = 

Tho stability or instability of the systen is demonstrated oy 
the use of a special plotting diagram in the complex plane surzeated 
by HW. Satche as a moans of handling aystems with time lag, and 
developed for this use by H. S. Taiem. This involves separating the 
trensfer function into two w rts. In the complex plane the first 
portion of the transfer function, the exponential variable containing 
the tins lag, plots as a unit circle as the complex variable p is 
made to take a contour enclosing the positive half of the replane. 

Af the remainine portion of the transfer function intersects this 
unit circle, the rocket motor can be unstable for large reduced tine 
lag; 1f it does not intersect the unit circle, the systom is 
generally stable, althoush the roots of the exponential coefficient 


in the positive half of the complex nlane mst be investigated. ‘this 
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Lattor rovatirencat econ te convenioentiy sceou dished by the aad of 
& Wyniuist o4aerran. 
The Gouuations for the feedbaek circuit are develowed end the 
oxidizer wul fuel trensier function reiuirsinents are determined. 
Two cesas of stable corbusticn end tyvo exses of unstovlios 
eomaistion are enalyzed. One unstsvie ames 13 3teSiliged oy tho 


ad.-ition of a feedbac). cirevutt. 
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SVBOLS AI) DIPTICIONS 


tine 
instantaneous value of the time lag 


pressure exponent of pressure dependence of ths 
processes taking place during the tine lace 


instantaneouns pressuze in the conbustion chanbor 


At = ot of the characteristic equation 
with the reduced tine as the independent variable 


pressure in the combustion chamber in steacy 
operation 


(p - p)/p =  froctional variation of pressures 
4n the combustion chamber 


instantaneous rate of injection, burning, and 
ejection of propellants 


common valus of the above rates of injecticn, 
burning and ejection in steady operation 


(fy - A) /n = fractional variation cf injection 
rate 


instentansous mass of sasss in the sconbustion 
shanber 


Me / m= gas residence time in steady cvsration 


O e +7 = total residence time of propoliants 
fn stoady onerstion 


absolute temperatures of commstion eases 
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jde 


av) 


Bis Fe 


if 


if 


if 


reduced aminlific. tion coefficient 
reduced anguinr Pra tency 
instantansous nixturo ratio 
mixture ratio in stexdy operation 


coefficient related to the variation of 
temerature with nixtuve ratio 


non-dinensional reservolr capecity for the 
oxidizer line 


non=dinaisivnal reservoir capacity for the 
fuel lines 


the devietion from unity of the mixture ratio : 


subscripts indicating injector end 
subdscrivt indicating ejector end 
subscript indicating oxidizer 
subscript indicating fuel 


t/O, = yeduced tine 


re 
ia. 


constant relating the mass flow rate to the 
pressure in the feed system Lines 


reduced tine Lear 


H] 


8/20 pA, O, =  {nertia parameter in the line 
p/2ayp = pressure drop paraneter 


ratio of functions in (a /dz) descrivine the 
response of the feedback circuits 


H+ KeKe”~ = a convenient croupin: of tems 





vid 


At, BS, ny a ne 


complex constants used in the trial 


solution 

(yp) =  gysten trensfor fumction 

8, (P)» g(r) = components of G(p) used in the Taten=Fatche 
Magram 

a, Co, D RFR, 5 = convenient aystem constants doerined in 


equations (16) through (22) 


Ly M = convenient grouping of systun const nts ond 
varieblos defined in eoustions (22) end (23) 


ne) = F, (p)L ~ Fe (pit = G3 /P + a > Gy cb Oy 


a convenient symbol for the feedbac!: circuit 
transfer function 





Ie TELINCTIS 


pocket motors often exhibit roucsh combustion cither as 
ehngring (low frequency oscillation below 100 cycles vex second) 
or scersaning (high frequaney oscillations severel hundreds of cycles 
per second). the first of these phenomena has bean investigated 
sufficiently that several suthors have developed equetions describing 
analytically the oscillation procsssss in rocket rotors ‘?), 10 30 
equations which are utilized in this paper are deduced ly L. Gcrocco 


under the following assum tions: 


1. ‘the processes of trensforméing the propellents into high 
torroerature gases require a cortain tins called the time les, T ; 
4% is the time after which the lionids nay be considered suddenly 
transformed into hot gases. ‘The rate of mixing of the provoliants 
denends essentially on the injection system used end recirculsticna 
effects, but probably not on chamber pressure. ‘The conditions, 
droplet size, temperature, vaporization and eventual chenical trans 
formation of the propellants after mixing depend considerably on 


chanber pressure. 


Ze ‘The time lag is considered the same for all particles of the 
propellants: the residence tine is considered the same for all 
particles, 1.¢-, @11 particles travel from the injectcr end where 


vurning oceurs to the exnanst end in the same rosidence tine interval. 





Fon) ye aC 2) a & & ff 7) unions, ie Sie em fon ae: 
=. Pho Lay cf ve Ppiation of tine Lac 30h Tvéssure 16 Bezumed 


J f (0) dt! = constant (1) 


$< 
whore p io the chamber pressure end f(p) 41s the tine rate of 


preperation for combustion. 


4 vrassuxs variations sre sufficiently slow throushout the 


co. wustion chembor that a uniform pressure can be used, p = 74). 


Se “Gemperature of the gases is aasumed to vary with time snd s--ce 


but not vith pressure oscillations. 


Ge Analysis of the stabilicine circuit itselr end its o3cilli- tions 
effecting the rocket notor are not investigated. 

Lie Oreos?” has denonatreated that a rocket combustion chanbar 
is intrinsicelly unstable, i.e., the buming phenonens itscll cen 
eenarate oscillations evon without variation in the injection rate. 
“lis oceurs because eny oscillations of the prossure will cause an 
osciliation in the tino leg by cquction (1). Thus if the tine lag is 
oscillatins, 1% is immediately seen that syhen the tine las is 
deecrensing, tze Dburnins of the particles that were injected later will 
cate uo ith the burnings of those injected earlier; this increasos 
the rate of Durning with rosvect to tho average. jJ4milariy, the 
Osits happens for en increasing tine lee, Thus, if the increased 


6P cocreased rate of burning should coincide vith a pressure 
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increase or decrease, self-excitod oscillations can exist resulting 
in unstable burning. 

For the bipropellant case assumption 2 that © g is the 
residence time of all particles means that each particle burns at 
the injector end and travels to the exhanst end carrying with it the 
temperature developed at the combustion Vast eis. Tims, the 
temmereture at the eaxhaust, Tope at the tims + corresponds to the 
nizture ratio of the propellants injected at the earlier tine 
t$-O,=t-T - 9 This can be put in teims of the reduced 
tine by dividing ty ©» giving: zs - S =1l., therefore, in effoct, 
we heave two reduced tine lags, S and1.0. I¢ is ths second time 
lag that introduces additional complications in the formulas for 
the bipropellant case by appearing as an exponential tem in go(p). 

The problem of stabiligzine a combustion chamber is simplified 
by the use of the TsieneSatche Diagram. If the system transfer 
function, G(p), is separated into a (p) = a? ? aa @5(p) = all 
renaining terms in G{p), then as the complex variable, p=A+ iw, 
49 made to take the contour enclosing the positive hal? of ths 
p=plane, stability can be conveniently determined if the G_(p) curve 
Lies completely clear of the unit circle, g,(p). A Nyquist Diagram 
of the denominator of g.(p) also is required to determine if there 
are any poles in ths positive half plane. 

For an umstable case in which g(r) does intersect the unit 


circle, @,(p), stability can be insured by designing a feedback 





a? 


ciroult transfer function which will cause the 65(p? curve to 

ove clear of the g,(p) curve, and still meet the requirenents 

of the ilyquist Diagram concerning poles in the positive heif plane. 
For a biprovellant rocket the teupcrature of the conbustion 

gases is a function of the local nixture ratio. As the mixture 

ratio operates away from the steady-stete velue, r, the fractional 


variation of the tempsratura varies as? 





m1, / r-r or, — te eG 
oe eh e K OF dr 


This defines K which turns out to be an interesting reramcter 
of the system. K may be aither positive or negative, although, 


as pointed out by L. Crocco, K is generally positive. 





Il. DIYTLOPMUU? OF BpaATIoNg 


Under the assumptions given above, the equation of the conbustion 
chamber for a bipropellant rocket was developed by L. Crocco: (See 


Meure 1 for schematic diagram of system). 


d : 
db + 2p, na] #9 RCMP = 


(F+A) Mel + (2-H) 4g]. + CA- #) (2) 


Following H. S$. @sien and considering the rropsllent to be fed 
by a centrifugal pur run at constant speed the following equations 


for the bivrorellant feed system are developed: 


SK 
Bea tot Clit a gis Lae Ms + Je Z (ut os) (3) 


ha 


Equations (2), (3) and (4) specify the transient behavior of 
the fuel and oxidizer sunply systems coupled with the combustion 
chamber when the reservoir capacity in the feed line, K, and Kp, 
are specified. The closed loop stabilization circuit can be 
conpleted by requiring Ky and Ke to be a function of P . This 


means that the reservoir capacity responds to nressures fron the 





esibustion chanher throng a suitable ampliffer end servo srstom. 


mus K, mi X- are assumed to satisfy a pair of tntogreeliffeorcntiol 


egus. tions, 
d : 
Fy (Gz) Ph =% is) 
i 
Fe (a> ) p — Ky (6) 


where Fy and 7, correspond to a ratio of functions (generally 


polynomials) in d/ds. 


Trying e. soluticn of tho form: io 
z 
gh a A‘e? Dye PP ae on Aiea 
/ Zz / PZ (7) 


where nowp= A +40 48 comlex, as are At, BY, Dt, L' and“. 
Uliminating the cormon factor 6°”, equations (2), (3), (4), (5) 
and (6) represent a solution of the system if the five homozeneous 


first degree equations beloy are similtaneously satisfied: 


E a“ apt 
ees fags oh-£.4 08 Bee - Bes Die o (8) 
Pat + [rs trghy) +a Be + | (Pot?) +d * |ut = 0 (9) 
6 ee Sal oP | & Mors /P op |b 


Peat + Jl + 2 (Pgt:) + Tgp |p [xCPgth)p + Jen? fe = 0 (10) 





ave 


(11) 


{ 
o 


B, (v) Alwj lex 


(12) 


{f 
fa 


Fe (p) AU we rit 


Since they are homogeneous equstions, they can be sinul tancousiy 
satisfied if, and only if, the foliowing determinant of the 


coeffictents of the above equations vanishes: 


Gn) t 
pt oe as ee Z 
ne g*) (q-4)¢ 


Fe [itz(Rord)p] © [tRed)p+ier] 0 


i o | [Kari 0 | [Gta pep 
E (p) © ° =| ° 
Fe) 6 0 : | 

where @ = Ht Mo Ke ? has been used for conveniencs. 


Mxpanding the determinant gives sauation (14) below; if p 
satisfies this equation, our assumed solution, equation (7), will 


be a possible solution of the systen. 





ow Ceo 


[pra-mgpie d(HreKpl['t <2) + Io pl] + 
oF rfirc(e+ efi tai +4)] + B(gtty[ita (e+ a) ] 4+ 
“Sn J [ite (B+ 4)] + 
-P,(g- 4 te (R44) ]§ + pe { fal a (lot @ 
2 -of 
nF [IteG+ 4) ] - PG) +R tH + PC mbH et 
peo rcp {asda PHIL ABATE # 
pret rep Gq erf[echt iI] + GLira(R+ s))} + 
pe? Hp (9-4) Cle 52) ie 
“ 4)[ a (Q+4 té(h +4] f Ft 
pe? Kup fl grt lars) IL/+ 2G tt 
se Ke) (gr 5) LAR AHI I + LLIB tBIt 


pe Pry (q+ 2)¢ Jo Ty) 2 Gp = © (/4) 





oo Sao 


Tor the systen to be stable, the roots p of equation (14) 
should not have real parts. The prodlem of stabilizing the rocket 
motor thus Decones one of designing ths closed loop tronsfor 
function, ¥,() * Fe(p). go thet this criterion 45 satisfied. This 
condition can ve investigated most easily by use of the Tsien-Satcke 
Magram and a yauist Diagram. 

If the complex variable p is made to take a contour as in 
Fisure 2 enclosing the positive half of the pepleane from - 00 to + CO 
end around a circle in the positive direction with infinits radius, 
each root of G(p) = O (aquetion 14) in the positive heif plane will 
cause G(p) to make a complete revolution about the origin. Soparating 
Gp) into two parts, & (>) end go(p), by dividing by the ccefficient 


; -Sp 
of the exponential term, e e &LVest 
GP) = 9,0) — Fe) 
_ ) [ber lAts pIL Bt Fo pF 
=e az OPT LE ee 
“HL I p+ (D+Seh p+ C+ReFI +EQL+ Kiem} 


where 


se = L+ x (Pp + 3) (16) 





B= L+t (e+ 4) (17) 
H 


Ds. Y 
Gan +(Pathe) Gt gn tc tie) + PePe) (E+ Eee + A it 


(i+ | ) | (18) 
Ds f (JetJg) (at) + (IgePotd Pe) (4 t2) + WePomyPe) (A+ 2K) (19) 


RE (Pp - Po) (e+ A) (20) 


S= «-K (Se Pood Pe) (22) 


x 


HS ce aK -»)/@- 1) (Bp + Sop”) + Je (Bp) + 


Je, v* | - Ke? [ a -1) p+ dp*) + Te(mp") % Ieop°]f 


(22) 


= { (H+ 2K + »)] @ ~ 1) (4p + gp”) + dy (4n") + SeJQr" + 


~ om? | (2 = 1) (ap + Sg”) + Fy (0?) + Spier? ] 
(23) 
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Note that in thsse definitions, A, 5, C, D, R, S, are constants 
of the systen, but that L and if include the variabdle (p). 

44th this arrangement of terms, G(p) becomes a vector with 
vortex on 8, (p) ani the starting point on g,(p). As p makes the 
stated contour of Figure 2, @; (p) decones a unit circle and go(r) 
prescribes various contours dovendent on its Oden: 

It 1s alse necessary to investigate the coefficient of ae: 
(this is ths denominator of @(>) as written in equation (15) to 
ascertain if 4t has any poles or seros in the positive half wilene. 

If the coefficient of 7 oP has 8 poles snd ¥f geros in the 
richt half peplano, then as p traces 2 clockwise contour as in 
Ficure 2, the coefficiont will trece res clockwise revolutions. 
therefore, for ataaility &(p) must meke r counterclockwise revolu= 
tions around the unit circle. Tims even though the primary condition 
of stability 1s that eo(p) lie completely clear of the unit circle, 
6, (»)s 4% rmet alsc satisfy the conditions dictated by the number 
of poles in the positive half plane of the coefficient of tno 


exponential torm, e ona 
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TIZ. STASILICY OF BiLHOPTLLAN? SPST 
WITHOUT F usDRBACK 
To investigate first the courled feed system and combustion 
chanber without the feedback circuit, set Kor) snd Ke(e) equal to 


gerde This reduces squoetion (15) to: 


(24) 


Gp = OOF jo 
~[nJo 5 $2 +(D +SeP)p+t Ct Re? J 


where the constants A, B, C, D, $ and B are defined in equations (16) 


throush (21). 


GtS8 ft: (Stable Gase, K= 0.2) 


i 
ge 
3 
me 


Leen = Os X = 1.0; F = 2.75; H = 0.233; T 
at, /de = 650; K = 0.23 Jo = 4.0 Jp = 0.86; P, = 0.53 
Pp = 1.0. 

Substituting these values into equetionsa (16) through (21) gives: 
A = 2.5; B = 2.0; G = 1.653 9 = 2.08; R = Gace} 


Gil? « 


tt 


S 


Thus, squstion (74) beconss: 


(p+ 0.9)(2.5 + a.8p)(2-0+4%0f) 


= (26 
- [0.32 p> + ().08 40.72 P)A C165 F9.2€ / 


Gye 





cr 


~Loe 


Now let p= 1) , where © is the reduced enguler freyuency. 
“austion (25) reduces to: 
6( da ( (3207+ 94400) 6 —IZ FEOF ES * 
= eee 
P ~J—0.32*+- (1.08 40.728 Jew +465 40.2677} 


Substituting velnes for © gives the Tsien-Satche IMegron for 
Go(p) plotted in Figure 3. Checking the denominator of g.(p) for 
poles gives the Nyquist Diasran rlotted in Figure 4. The nunsrstor 
is plotted in Figure 5. 

Fron the figures 1t is evident thet this case is stable since 
the g.(p) curve is completely clear of the unit circle end mekes no 
loops about it. Also, the denominator in the Nyquist Diesran is 
proven not to have any zeros in the positive helf peplene, since it 


does not make eny revolutions about the origin. 


CASS II: (Stable Case, K = 0.3). 

Letn = 0.1; % = 1.0; B= 2.75; H = 0.233; 2, = 4500°% : 
a2,,/ar = 1000; K = .3, Jo = 4.0; Jp, = 0.83 2, = 0.5; 
ve ~— 1.0. 


Substituting these values into equations (16) through (21) givass 
Bee 8 2.5: B= 32.0: GC = 1.650: D = 0.26; BMS Ooo Peer 
Letting p = 41) and substituting these values into equation (24) 


gives: 





be -32w + (544-6), — 4A 4pore 45 
G@=e - aia Ce (27) 
[0.32 d*+ (0.36 those" Jew + US 4 038°) 


Azein, substituting values for © gives the Tsien=Satche 
Diegram for &(p) plotted in Figure 63 the liyquist Marran of the 
denominator is plotted in Figure 7. Ts nunerator is the acme as in 
tas0 I (Meure 5). 

Chaneine K in this second case has rasulted in a radicel chime 
in the TsieneSatche Miseram. Ths curve &5(p) now makes two counter 
clockwise loops about the unit circle, althouch 1% does not intersect 
the curve @(p). Nevertheless, this is still a stable case since 
the denominator of &3(v) mass two clocikvise revolutions about the 


Origin in the Tyquist Diesran. 


CASS Lil; (Unstable Case, K = 0.2) 


Let mn = 0.5; A = 1.0; F = 2.75; H = 0,233; 2, = 4500°K; 
afg/dr = 650; K = 0.2; Jy = 4.03 Jp = 0.83 PL = 0.5; 
ie a 1.0. 


Substituting these values into equations (16) through (21) gives: 
om 2.6: BR = 2.03 9 «= 3&.65; D = 5.72; (Rm 0.2} 


S = 0.72. 


Letting p = 34) and substituting these values inte orution (2-4) 


“givess 





Gl dw (-3.20° + 10.80 dé — 73-20" + 2:5 (28) 
= C8 ee id 
“ ~[-106 2+ (572 F0.72E7P) CO + (BF FO.2E°°)] 


Substituting values for “~ end plotting gives the “sioneSatche 
DMiagran of Figure 8 for @o(p). Tho behavior of the denominator of 
@5(p) is shown in the plot of Mgure 9. The numerator is vlotted in 
Tigzore 10, 

Inspecting the figures it is seen that the B_(p) curve now 
intersects the unit circle and thus the rocket motor cam be unstable 
for large reduced time lags. This shift to instability is prinorily 
the result of the change in the value of n from 0.1 to 0.5. ‘rus, 
tn", the pressure exponent of prossure dopendence of the proceosses 
takine place during the time lag, T , is shown to be an imortent 


systen parencter. 


CASS (Unstable Case, K = 0.3) 
etn = 0.6; % = 1.0; Fe 3.76; H = 0.238; Sh = abe0 7; 
an / ar = 10003; KE = 0.35 dg = 4.03 Jpg = 0.8; PL = 0.53 


ie a 2.0. 
Substituting these values into eyuations (16) through (21) gives: 
mee cong B= 32.03 C€ = 3.80: Dae 50; F = Clas 


5 = 1.08. 





~1o< 


Letting p = 4 Wand sudstitutine these values into 


equation (24) cives: 


_ dw (-3.2°+ lo- we -—/22m7 4 25 


a oe eee 
Gp-e -[- 16 W*+ CES #408 a“ Jew + (3-540.3e°%] 


Substituting values for © and plotting the values of Bo () 
gives the SsleneSatche Diagrem of Fieure ll. The Myauist Macran of 
the denominator is plotted in Fieure 12 and the numerator, which is 
the sane as in Case III, is plotted in Fievre 10. 

in this case chancineg K fron 0.2 to 0.3 did not have the 
spectaculsar results as in the stoble ease, although the g(r) curve 


does have a more commlicated shape. 
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IVY. STABILIZATION OF A BIPROPELLANT SYSTS: 
BY A FeaD ack Cincul® 

Taurine to the problem of stabilizing an unstable burning 
rocket motor, it is obvious that a feedback transfer function, 
By (p dL + Be(p)M, 48 required that will move the eo(p) curve 
completely clear of the unit circle for smali Siives of DD) plus 
satisfying the requirements of the Hyquist Diagram for all values of 
uw). Since the B5(p) curve can only intersect the unit circle in 
the region of am-ll values of i e the ae terns in the denorie 
nator of Ep) — (equation 15) — can be expanded in a power series 
of (cw) and higher order terns neglected. Also, define a feed 


system transfer function, ¥(odk + Fo(p)H, as a series in povers of p: 


F(pIL + Fy(p)M= ~ (p) = ay/p t a +t apt aon” (30) 


Substituting these two chanzes into G(>), equation (15), and 


letting p= 40 gives: 


dws /~ cw ea (Cw 
Gw-¢ Z [cw t0-n) [Atco [ B+ Je ( ee me 


{rE i. Gur"+ +X] - co + + JJew +C#eRG- (w+ rT, is 


Mow clear tho denominator of go(p) of imeginary teras by 
multiplying both numorator and denominator by the conjugate of the 


denominator giving: 





~li&e 


d (cw) = - [ban yjln Jo e+ = -S+ G2] -[D-R+s + 2, }[3.A(-n) +Js B(r-n) + BA]} u)* 
2 jz[n F348 -st+aj[ct R+4o]-[D-RrSta;]* (éa)*4 [c+Rt a] 


SILA +5 B+JeJo(-n CIR + aor” + [C+R4+4,][BAC-7) ] 
f2[nby +&_S+a, ][CtR+a.J—[D -R+S4G,) (ca) + [C+ R+a,1* 


LEE AU 1) + TeBG-n) + BAILOR +40] - DR +54 aL BAC © U0, 
jal ye 45 +8 -S+ 4, [1G +4.]- [O-R#34%] Jes)” + [CHR A. | 


This is the approximated equ-tion for the go(p) term of the 
system trensfoer function and holds for small ™ only. Terms of 
order (4 ~ ‘7 and higher order were neglected in this develoynent. 

Now it 18 necessary to specify the shape of the stabilized curve, 
G5(p), near the unit circle and dorive the values of Ay. 22 S a5 


that will meet these roquirements. Let this stabilised curve ve: 
‘ 2 3 a 
Go(p) = V+ (1H) + TMH)” + 24) (33) 


Aquate this to &-(p) in equetion (31) and clear by multipiying 
the denominator by equation (33), transposing the numerator terns 
to the left hand side and collecting even and odd powers of (iw ) 


togethers 





= Ou: 


fw[n LJ, + $-s]+xD-R+5] + ¥ [c+ R]+ Wart Xa, + Yao t Za4§(u)"+ 


[LA +5 B+I,3.(1-n}] Cw)? + JW[CHR]+ ABG-n) t+Wa, +X a4 f=20 (ca) 


iw -8y +x [ni s,+8-s] + VID-R+5)+Z [e+ 1+ Xa, +Ya,+Z ao} (Cu)? 


[3.1 cu? + 2 w[p-R+5I +X [CR] + [LAG-n) +Jp BO-) +BA] + 
a (35) 


Wa, #Xa,t YG, {wy + We, = © 5) 


Hote that in this development terms of order (3w)* and higher 
were neglected. The coefficient of exch (4W) term mist be indevend= 
ently zero for thess e.untions to be satisfied. Thorefora, setting 
the coofficient of esch term, 1.6., (LW en (QW), (hw ia 


2 
(42 )” to zero end solving for the mknowe, Bipe Boe Ae By Gi vest 


Waly = 0 
Hag + Xa == (len) - wic+ 24 


Wa, + Xay + Ya. = - | g,A(t=n) + Seb (len +34] =‘! [ mR+ 3] -x[o +R] 


Wag + Xay + Ta, + Jay == [IgA + Jeb gly (1en)| ~u [niga + " -3] ” 


j 


~z[n-nes] -r[ce¢ 


ad 





=20— 


fos3s four equctions can be solved simultaneously to sive 
explicit relations for 2.4, Soe 15 &, in terms of system const-nts, 
ont thas 3, (po) Bp(o Mit, the feodbect: system trensfer function, is 


datorained, 


Solving for the unlmown 8.4. ao» %» a> gives: 


2.1 = 0 


a= - Loren - (C+ 8) 


a = -] 5,4(2-n) + JpB(l-n) + Bal - (DRS) + (C+ BR) (Xx -% )+ 


[ ap(ien) | 


aye [ JAGlon) + JgB(l-n) + nila ee | ni gJ¢ As - s| 


[ mao) -X5] -[o+ 2H - 22] 
(37) 


For amll © the approximated g.(p) equation (32) can be used 
with the feedback transfer function equation (30) to plot the 
@sien=-Setche Diagram. For large “~ the general equation (15) mst 
oe used in plotting this dlagrams also the Uycuist Diagram of the 
donominator of B5(v) nust be investigated. 

To Jllustrate the proccdure outlined above, consider the unstable 
ease discussed previously where K= 0.2. Let the required stabilized 


curve of g@5(p) be doscribed as follows: 





J 
me, 


‘ ‘ any 
G_(p) = Wt X(1w) + Yio) + ato) = - 24 6 +819. (38) 


This arbitsary selection of the g,(p) curve 1s made to plore 
the stabliized curve in the aroroxzimete position of the know stable 
curves of Cases I send II. ‘thos, Yee 23; X= - G 
Also, A= 2.53; B= 2.03 0 = 3.653 De 5.723 Re 0.03 $= 0.72 (other 
values as Sase Iii). 


Substituting these values into equation (37) gives: a, = 0; 


Bg == 263 & = - 5.64; a5 == 1.06. ‘ith these values, 
E = = 2.6 + 1.060" = 5.64 1. 


Using equation (32) and the above values the following expression 


43 obtained for Bp): 





13 - (5.391 O~ + C(/2.46W) " 
§2P) a. a as 2 2 
4/53 — 3.28 W 


Solving this equation for various values of “ gives ths bslow 


$ablo of results: 








zs d2(P) We Xu) (ced? + Ze? 
O Or —2.00 

Es 198 - ¢ 03 Ee a 5: 

A —-($6- 62.25 -(5 ~< 2.G0 

3 -/22 - ¢ 3.55 -/.12 - CLEG 

z £0.97 —-¢ &.7 o—-¢ Zo 





ee. 


This table shows that as © becomes groater then 1/3, the 
values of g,(p) have too great an error. However, from Mguro 13 
4% 48 apparent that the curve hes bean atabilissd satisfactorily 
for small y). 

It is now necessary to investigate the dehavior of Bo(>) 
using the exact equation (15) which includes tie feedbeck circuits 
substituting the values determined above into equetion (15) rosults 


ins 


Sp | -/3.20742.5 4+ ¢ (3.2wW? H/0.Pw 
Gp=<« P : ay : (40) 
-[-o.s¢w + (0.12 40.72E“ Jcus + 105 40.2E°% 





Substituting values for “© and plotting &_(P) gives ths 
Ysien=Satche Diagram of Figure 14. ‘The liyquist Diagren of the 
denominator is plotted in Fieure 15; the numerator is the seme as 
Case III and is plotted in Meure 10. 

Pron Figure 13 4t 4s evident that the stabilized curve is 
completely clear of the unit circle for small ~ . However, for 
large © the addition of the feedback circuit has caused the g5(p) 
curve to chance in the Tsten-Satche Diacran from no loors about the 
unit circle in Case TII to one with two counterclockwise revolutions 


adout 8, (p) in the stabilized curve of Figure 14. WNeverthcolessa, 





ma 


she system is st411 st-bvie since tho tyguist Dlevran shoys that 
the denominator of &.(p) males tivo clockwise revolutions about the 
Origin. 

Thnas an unstable burning rocket motor has been stabilized by 
en aporopriate closed loop feedback system. ‘The cireult chosen 
to stabilize Case TII was a combination of three terms of a 
differentiating circuit and one tem of an integrating circuit. The 
coefficient of the integrating circuit turned out to be zsro 
indfeating this tyne circuit need not be used. 

In order to show the nature of the feedback function, b (>), 


set Fp(p) = 0 and solve for F, (>) fn equation (30). 
2 ; 
(>) = (a_, /pbay+a, ptop”) eS Ay /p +4, + AyD + A.D (41) 


nvalinuating L by substituting the values of constents used 
72 


previously, eoproxinating ei? by a pover series in (12 ) and 


neglecting hisher order terms gives: 
a) 
B= —- 10.59 (4@ )” = 2.33 (14) (49) 


iquation (41) can nov be evaluated by using the known values of all 


terms, resulting ine 


Pi(p) = fe 1.06(4@)” ~ 5.8444) - 2.601/ [+ 10.59(4) + 2.53(4 )] 


2 
== 1.12/ (14@) + 2.57 = 12(42) + 54.77(4~) (43) 


4 ’ 


muerefore, to our order of approzimetions 





= + 3.57 
*e (At) 
Ay se 12 
Ap a <p> 54.77 


It is now nocessary to choose 2 relationship for ¥, thas will 
encroach zoro for larze @ and still retain the desirad stabilizine 
characteristics for smell WwW . Thus the following expression, which 
is arbitrarily chosen, describes a circuit that wili stabilise the 


eonbustion chambers 
2 
7,(p) = (-1.12/p) (1+ Bp / 1 + Bop” + Bap”) (45) 


fo evaluate Bye By» ae in tarms of the constants in equction (44), 


aouste equation (45) to (43) givine: 


By = — Bend 
B, = - 10.7 (46) 
Bs = + 23.9 


Thus, 


Fi(p) = (-1.12/p) (1 - 2.29p/ 1 - 10.7ip® + 23.9p") (47) 


This ecuetion for the closed loop circuit gives the desired 
stabdilizetion to the rocket notor for low frequency pressure 
variations in the buming process which could have caused combustion 


instability. 
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